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Abstract 

A semi-analytic method to compute the first coefficients of the renormalization 
group functions on a random lattice is introduced. It is used to show that the two- 
dimensional O(N) non-linear <r-model regularized on a random lattice has the correct 
continuum limit. A degree k, of "randomness" in the lattice is introduced and an estimate 
of the ratio Arandom/ ^-regular for two rather opposite values of k, in the <r-model is also 
given. This ratio turns out to depend on k. 
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1. Introduction 

The lattice has extensively been used in numerical simulations and has proven to be 
a unique tool for studying the nonperturbative aspects of a field theory [1,2]. However, 
some of the Monte Carlo simulations have provided only a qualitative explanation or 
at most the order of magnitude of the quantity under study. This problem holds for 
example for the work done to extract expectation values of composite operators (topo- 
logical susceptibility, gluon condensate, etc.). In these simulations, the physical signal 
is hidden by mixings whose perturbative expansions are badly known since few terms 
have been calculated (and the series are at best asymptotic). The random lattice [3] 
could be useful in this respect. Indeed, the scaling window on random lattice QCD is 
strongly shifted towards small values of the coupling constant [4,5]. This fact may allow 
a better control on the perturbative expansions which mask the nonperturbative signal 
(perturbative tails, mixing with other operators, etc.) as well as on the non-universal 
terms in the scaling function. Obviously, a careful analysis must be performed in each 
case. In particular, for QCD, the onset of asymptotic scaling should be deeply studied 
on a random lattice. 

The coefficients in the perturbative expansions should not be larger on random lat- 
tices than on regular lattices, otherwise the advantage of working with smaller coupling 
constants is lost. One of these coefficients, the first term of the perturbative tail of the 
plaquette operator, can be easily computed by using the equipartition theorem. On the 
4 dimensional random lattice, it is 0.1364(iV 2 - 1) for SU(N) QCD [6], while on the 
regular lattice it is 0.25(iV 2 - 1). This is good news. However, an explicit calculation in 
each case must be done in order to decide whether the expansion has a better behaviour 
or not. 

A question raises naturally at this point: how do we compute the coefficients of the 
expansions? On regular lattices, the perturbative expansions beyond 4 or 5 loops become 
rather awkward due to the fact that vertices develop a huge quantity of rotationally non- 
invariant terms [7]. If one managed to get rid of these terms, thus having vertices with 
the same size as in the continuum, computations would become much more feasible. 
The almost-rotational invariance of the random lattice could be exploited in this sense 
[8]. Indeed, in contrast with the regular lattice, the random lattice displays a higher 
degree of rotational invariance. In particular, from the nonperturbative point of view, 
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numerical studies of nonabelian gauge theories have revealed the absence of a roughening 
transition [6,9] from strong to weak coupling regimes, thus opening the possibility for a 
wider scaling window. 

Despite these potential advantages, the random lattice has been barely exploited in 
numerical simulations. Especially, some of its properties have been little checked. For 
instance, the continuum limit of models defined on random lattices has not yet been 
completely clarified. For statistical models with power law scaling behaviour for the cor- 
relation length, the Harris's criterion [10] predicts whether universality on a lattice with 
impurities holds or not if the specific heat critical exponent a does not vanish. For the 
two-dimensional Ising model (an example of vanishing a) recent simulations [11,12] have 
shown that the critical behaviour depends on how one introduces the "randomness". 

Asymptotically free field theories regularized on a lattice provide another important 
kind of models. In this case, the correlation length shows an exponential dependence 
on the inverse of the coupling constant. An extension of the Harris's argument [13] 
supports the scenario where the regular and random lattice regularizations share the 
same continuum limit. However, this is rather a physical consistency argument and an 
independent check is worthwhile. 

In this work, we give an analytical proof of the universality for a field theory re- 
gularized on a random lattice. We have performed a weak coupling calculation of the 
first coefficients of the renormalization group functions for the two-dimensional O(N) 
non-linear <r-model regularized on a random lattice. We have followed a semi- analytical 
procedure which will become clear along with the explicit calculation explained in sec- 
tions 4 and 5. The action for the O(N) non-linear <r-model in the continuum is 

S = j-Jd 2 x(dJ) 2 . (1) 

— * 

The arrow on the field </> stands for its N components which are constrained by the 
condition (j)(x) 2 = 1 at every space-time point x. In Eq. (1), g is the bare coupling 
constant. We have chosen this model because it has some properties in common with 
QCD. It has a rich topological structure for N = 3 [14], it presents spontaneous mass 
generation and is asymptotically free [15,16,17] in two dimensions for N > 3. This last 
property is the one we are interested in. 

We have defined the random lattice as in reference [3]. We have also included 
a parameter k, to describe the degree of "randomness" of the random lattice. This 
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parameter is essentially the minimum distance between sites. We have checked that the 
same continuum limit is reached for several lattices defined with different values of k. 

The previous computations also provide the ratio between the renormalization 
group invariant A parameters Arandom/ ^-regular- This ratio is an essential ingredient if 
asymptotic freedom is to be checked. We have shown that this ratio depends on k. This 
is an important result, mainly if, as originally proposed by T. D. Lee and collaborators 
and afterwards used in simulations with fermions, one is interested in averaging among 
different random lattices. Indeed, this k, dependence means that one has to average 
only among lattices with the same k. One can also take advantage of this dependence 
to improve the asymptotic scaling of the theory in the way we will discuss in section 6. 

Before ending this introduction, we make some comments on the feasibility of a 
simulation on random lattices. Any numerical study of a theory on a random lattice is 
necessarily more involved than on regular lattices. Firstly, one has to store the infor- 
mation about links, plaquettes, etc. which can be rather memory- consuming. Another 
problem is the parallelization of the updating algorithm. If one can use a cluster algo- 
rithm, then the implementation is trivial and only a quite small increase in CPU-time 
per updating is seen (due to the larger coordination number on a random lattice in 
any dimensions). In QCD, instead, one usually makes use of local updating algorithms 
that need be parallelized to become efficient. On a regular lattice, this parallelization is 
well-known and can be implemented by gathering all sites in two sets, called colours. On 
a random lattice, a similar task can be done [18]. The number of colours on a random 
lattice is larger than on a regular lattice, therefore an unavoidable slowing down in the 
performance is detected. However, this slowing down seems not so drastic to prevent 
the use of random lattices. 

The paper is organized as follows. In section 2 the random lattice, the parameter 
k, and the action for the <r-model are introduced. This action on the random lattice is 
written in a convenient form for the subsequent calculation in section 3. In section 4 
the general expresion for the bare 1-loop propagator on the random lattice is written. 
In section 5 we perform the numerical calculation of this 1-loop expression and of the 
corresponding renormalization constants. From them, the first coefficients of the (3 and 
7 functions are computed as well as the ratio between A parameters. The discussion of 
these results and conclusions are presented in section 6. 
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2. The random lattice and the action of the cr-model 

In order to construct the random lattice, we perform the following two steps. Firstly 
we fill the volume V of the lattice with Af sites and then a triangularization process [3] 
is applied to these sites. 

To fill the lattice volume, we defined a minimal distance between sites by demanding 
that there are no two sites closer than y/V/Af / k, where k is some parameter to be fixed. 
The larger the value of k, is, the more irregular the pattern of sites looks. Instead, for 
small values of k, the sites tend to be more evenly arranged and to stay far from each 
other. Hence, such random lattices with small k, can be compared with the lattices 
constructed in reference [19] by using the eigenvalues of complex random matrices. 
Moreover, the parameters defining the structure of the random lattice (i.e.: link lengths, 
plaquette areas, etc.) show less variance for small k. 

The smaller possible value of k, is 1; at this value we recover the distance between 
sites on the regular square lattice. 

Once the sites have been placed on the lattice, we proceed to construct the trian- 
gularization [3]. We followed the method of reference [9]. It consists in joining sets 
of three sites to form a triangle with the only condition that the circle circumscribed 
by these three points does not contain any other site. The sides of that triangle are 
the links joining the three sites and the triangle itself is a plaquette. This construction 
fills the whole lattice with no overlapping among the triangles. We imposed periodic 
boundary conditions. Thus, to construct the triangles and links, we considered the two- 
dimensional lattice as a torus. In figure 1 we show two Af = 100 sites random lattices, 
with k, = 1.3 and k, = 100. 

One can define also the dual lattice. Its dual sites are the centers of the above- 
mentioned circles. It is clear that any link is the common side of two triangles. Thus, 
every link of the random lattice must be surrounded by two dual sites. The line joining 
these two dual sites is called the dual link. Therefore, every link is associated with a 
dual link. 

As soon as the lattice has been constructed, one can devise tests to check the 
triangularization. The first and easiest one is to verify that the number of triangles 
(links) is equal to twice (3 times) the number of sites [3]. Other good tests are the 
integral properties [20] 
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£vS = o, Evs^ = 2^. (2) 



In this equation, is the link vector joining the sites i and j and the matrix \{j is 
defined as follows: 



f 0, if z and j are not linked; . . 

lJ \ s ijl^iji otherwise, ^ ' 

where l{j is the length of the link vector and S{j is the length of the associate dual 
link. A third test consists in demanding that the action has only one zero mode (see 
next section). This last condition is very exacting and if some links lack, then spurious 
extra zero modes can develop. 

After all of these preliminaries, we can write the action for the non-linear <r-model 
on a random lattice as [20] 

5£ = iE A «(£-&) 2 - ( 4 ) 

9 i,j 

g is the coupling constant. i,j, ... denote sites and <f>i means the value of the field <f> 
at the site i. The arrow on the field <f> stands for its N components which, as in the 
continuum (see Eq. (1)), are constrained by the condition <f> 2 = 1 at every site i. For a 
regular square lattice, \{j is 1 for linked sites and zero otherwise. Hence, Eq.(4) becomes 
the standard action when it is considered on a regular lattice. One can show that the 
naive continuum limit of Eq. (4) is the correct action for the model in the continuum, 
Eq. (1). 



3. The Feynman rules on a random lattice 

In this section, we will take advantage of the constraint (j) 2 = 1 to rewrite the action 
in a more tractable form to be used in perturbation theory. 

By using the constraint, the field can be written as (j) = (7?, y/l — 7r 2 ) where 7? is 
a (N — l)-component field. The single components of the 7? field will be denoted by 
capital letters A,B,... as in ir A . In terms of 7?, the action (4) can be written as 
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/ 13 (5) 
where some constant terms have been left out. The last term comes from the jacobian 

— * 

of the change of variable from X></> to X>7r in the functional integration. The matrix Aij 
is equal to 

= A^, if i ^ j; 

A« = -£,-V < 6 > 

Notice that ^ • Aij = for all i. 

We add an interaction with an external field h to regularize the IR divergences 
which will appear. This interaction term can be chosen as 



where a is the lattice spacing on the random lattice defined as a = y/V/Af and a 
constant has been omitted in the expansion. 

Eqs. (5) plus (7) is the action we will use. It is convenient to diagonalize the matrix 
Aij which is equivalent to diagonalize the quadratic sector of the action. Let us call 
— £ a the eigenvalues of Aij and B{ a the relevant orthonormal matrix of the change of 
basis. It can be proven [21] that for the matrix Aij, all — £ a are negative except for one 
of them which vanishes. We have noticed that the fact that there is only one vanishing 
eigenvalue is rather sensitive to the correctness of the triangularization process, (see 
previous section). 

The index i in the matrix of the change of basis B{ a stands for the space-time 
position (sites), while the index a stands for the diagonalized basis position. Notice that 
this procedure is like a Fourier transformation and the diagonal basis is equivalent to the 
momentum space. Therefore, the B{ a matrix is equivalent to the Fourier transformation 
integrand exp(za;ip a )/(27r). Indices in the diagonal basis will be denoted by the letters 
a, 6, . . . as in £ a or B{ a . Although we use the same letter, a, to indicate the lattice 
spacing and a subindex in the diagonal basis, there is no possible confusion. The action 
(Eqs. (5) and (7)) written in this basis (in momentum space) is 



S L +S h =±- Y.Ua + ha^l (8.a) 

a 

(«•») 

a 

+ ^-{ \] ha 2 (TT a ■ 7T 6 )(7T C • Tf d )Bi a BibB ic Bid 

- ^ Aij(7T a • 7T 6 )(7T C • Tf d )B ia Bi b Bj c Bjd} (8.c) 

l,.7,<l,&,C,<i 

+ 0(tt 6 ). 

Now, one can easily extract the Feynman rules from this lagrangian. Eq. (8. a) gives 
the propagator 

8 AB a 

A AB = — 9 — (9) 

while Eq. (8.6) is a measure term, similar to a counterterm, with an extra power of the 
coupling constant g. The four-legged vertices are in Eq. (8.c). One usually performs 
the sum on i and j to get the Dirac delta indicating momentum conservation. Here, 
this step cannot be done. Vertices with more than four legs will not be needed in our 
computation. 



4. The one loop propagator 

(2) 

In this section we will compute the Green function at one loop. The tree 

order is proportional to the inverse of Eq. (9) 

r£e order = kh+h). (10) 

y a 

At one loop we have to consider the two Feynman diagrams of figure 2. The first 
one, figure 2. a, comes from Eq. (8.6). It is the measure and gives — 1/a 2 . Eq. (8.c) 
contributes with the tadpole shown in figure 2.b. This tadpole yields three terms, Ti, T2 
and T3. The first one is proportional to h: 

T\ = ^(^^ + 1) J2 B ^ B ^ A c B - (11) 

i,c 
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Regularizations usually preserve the momentum conservation for every value of the 
cutoff. In this case, Eq. (11) turns out to be proportional to 6 a f,. This fact hap- 
pens for example on a regular square lattice. On this lattice, the expresion A AB = 
^2 c Bi c A AB Bi c is independent of i and the subsequent sum on i in Eq. (11) provides 
the S a t, because of the orthogonality of the matrix B{ a . On a random lattice, we can 
assume that this equality is approximately true and rewrite Eq. (11) as 

Tl = h (^Y^ + 1 ) SABA oS ab , (12) 

where 8 AB Ao is the z-independent value of A AB . Now, we need a simple definition 
for this Ao- To this purpose, let us remember that among the eigenvalues £ a there 
was a zero mode. This is due to the fact that ^ • Kij = for every i. Therefore, 
the eigenvector associated to this eigenvalue must be of the form (l/v^AT, 1/a/AT, . . . ). 
Then, using the previous definition of A AB , we can define Ao as 

The validity of this assumption must be explicitly verified. For a random lattice of 40 2 
sites, Eqs. (11) and (12) differ by less than one per mille when using Eq. (13). 

The terms T2 and T3 do not display a linear dependence on h. The first one is 
proportional to 

T 2 oc^A^A^A. (14) 
Under the assumption of j-independence for A AB , this term can be approximated by 

T 2 oc S AB A B i<* B ib S A *i' ( 15 ) 

i 3 

and this last sum over j vanishes. Hence, without the assumption, one expects that T2 
is negligible. Indeed, an exact calculation of the Eq. (14) on a random lattice of 40 2 
sites, shows that it is less than l°/oo of T3. Thus, we will not consider T2 any more. 
The T3 term is 

T 3 = -^Y, Ai i A c B BiaB jb B ic B jc . (16) 
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This equation needs be written in a more manageable form. On a regular square lattice, 
Af B = BicA^ B Bj c is constant if i and j are linked sites. On a random lattice we 
can assume a similar behaviour and call 8 AB Ai such constant. If we further assume that 
An can also be substituted by an z-independent constant, —A (on our random lattice 
this is true within a 1%; on a regular lattice A is exactly equal to 4), then Eq. (16) can 
be easily rewritten as 

T 3 = ^S^Sa^A.Ca + A(A - AO). (17) 

Below, we will give explicit expressions for Ai and A. Now, let us prove Eq. (17). 
Splitting the sum in Eq. (16) in a term for i ^ j and another for i = j and using the 
previous definitions for Afj B and A^ B , we get 



T 3 = --S^^^A^iaBjb + Y, A "( Ao - A i)*iA + S A^uBiaBa). (18) 

i^j i i 

Now, recalling that Ai and Ao are constant, we obtain from the first and last sums 
AijBiaBjb = -£ a t>ab- The second term in Eq. (18) just gives A;;<5> a6 (Ao - Ai). 

Collecting all of theses pieces, we obtain Eq. (17). 

Let us write explicit definitions for both Ai and A. Recall that the eigenvector 

associated to the zero eigenvalue is (1/ y/N , 1 / y/N , . . .). By using this expression in Eq. 

(16), we can rewrite the £ a -independent term in Eq. (17) as 

^ AB (1-WA ). (19) 
Comparing this equation with Eq. (17) yields our definition of Ai, 

. . ha 2 . 1 . 
Aj = A (l + -=-)- (20) 

The simplest definition for A is A = — ^^Aii/A/". However, Eq. (20) suggests this 
second definition: 1/A = —{Y2i^/A.a)/Af. Let us call Ai and A2 these two possibilities. 
We will choose the definition which makes Eq. (17) closer to the true value, Eq. (16). 
As we will see in next section, we will work with two values of k: k, = 1.3 and k, = 100. 
For k, = 1.3, both Ai and A2 coincide within errors; but for k, = 100 they yield different 
numbers. Using any of them, the agreement between Eq. (16) and (17) is quite good: 
within ~ 5%, but A2 shows the better agreement. Hence, we chose A2 as our definition. 
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Moreover, an explicit numerical study on a random lattice of 40 2 sites shows that the 
value of Eq. (16) for a ^ b divided by the corresponding value for a = b is ~ 10 -4 . This 
favours the presence of the 8 ab in Eq. (17). 

(2) 

Collecting all pieces, the 1-loop contribution to is 

lf\ oop = SAB ^(-^ + K^ 1 + l)Ao + + 1(1 - k 2 A„)) . (21) 

The procedure to know the numerical value of Eq. (21) will be the following. We 
will construct an explicit random lattice and calculate the matrix A{j . The eigenvalues 
of this matrix will be used to compute Ao . The diagonal of A{j will provide A and the 
value of Ai through Eq. (20). This will be done in next section. 

However, there is one more problem to be solved. To compute the renormalization 

(2) 

constants, we are interested in the divergent terms of V) . If the behaviour of the 

' o bare 

eigenvalues £ a near the continuum limit can be written as a power series in p 2 , where 
p is the momentum in the continuum, one can reasonably assume [22] that Ao as a 
function of ha 2 , is expressed as a series: 

A (ha 2 ) = ^{ha 2 ) n c n + \n(ha 2 ) ^{ha 2 ) n c' n . (22) 

n>0 n>0 

To compute the renormalization constants we are interested only in the Co and c' terms. 
To extract them we will calculate Ao for a given random lattice and for several h. Then 
we fit Eq. (22) to these values of Ao(ha 2 ). Obviously, we need the values of Ao(ha 2 ) 
for infinite size lattices. We observed that the results obtained with a lattice of 40 2 
sites were essentially the same than with larger lattices if ha 2 > 0.1. That means that 
our finite size parameter hV is of the order of ~ 160 or bigger. Moreover, to make the 
fit reliable with only the first terms of the series in Eq. (22), we cannot consider large 
values of h. As a consequence, we are forced to work in the window ha 2 £ [0.1, 0.9] when 
doing the fits. We used least-squares fits. Alternatively, one could also make use of a 
linear system of equations: the values of Ao(ha 2 ) for m different ha 2 to determine the 
first m coefficients in Eq. (22). We obtained similar results than with the least-squares 
fit. 
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5. The numerical results 



In this section we will explicitly construct the random lattice and make use of the 
definitions of Ao and Ai in Eq. (13) and (20) to know the numerical value of Eq. (21). 

We used random lattices of 40 2 sites. This lattice size is enough to obtain rea- 
sonable results as discussed in the previous section. We chose two different k, para- 
meters: k, = 1.3 and k, = 100. Then, for each value of this parameter, we construc- 
ted 40 random lattices and determined their respective matrix A{j. These matrices 
were diagonalized to compute Ao by means of Eq. (13) for the following 16 values 
of ha 2 : {0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85}. The 
value of A was also extracted. The diagonalization of the matrix is performed with 
a NAG subroutine and it takes ~ 23 CPU seconds for our lattice size with a CRAY 
C-94/2128. 

At this point, we have 40 independent determinations of Ao for each value of ha 2 . 
Then, the final result for Ao(ha 2 ), to be used in the fit, is extracted simply by computing 
the mean and statistical error from this set of 40 determinations. We noticed that there 
is no correlation among these 40 values, so we assumed a normal distribution to compute 
the errors. 

These 16 values of Ao(ha 2 ) both for k = 1.3 and k = 100 are shown in table 1. 
From this table it is clear that the regularization depends only on k. This statement is 
supported by the fact that the determination of Ao is rather insensitive to the particular 
random lattice used, as the errors indicate. 

Now, the 16 values of Ao(ha 2 ) are fitted to the first coefficients in the series of Eq. 
(22). To accept the result of the fit, we demanded that the values obtained for these 
coefficients remain stable when we remove or add some of the 16 fitted points. The 
value of the % 2 is not a good test of reliability for the fit because we are not fitting 
measured quantities with a gaussian distribution of errors. Indeed, we always obtained 
a very small % 2 . The best result was obtained when fitting five coefficients in Eq. (22). 
Here we report all of these five coefficients, but recall that only the first two (co and c' ) 
are interesting for us: 
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A {ha 2 ) = 0.2825(29) - 0.07964(83) ln(ha 2 ) - 0.0255(30) ha 2 

+ 0.0141(35) ha 2 ln(ha 2 ) - 0.0007(20) (ha 2 ) 2 ln(ha 2 ) + 0((ha 2 ) 2 ) (23) 
A = 3.887(7), 

for k, = 1.3 and 

A {ha 2 ) = 0.246(11) - 0.0822(30) ln(ha 2 ) - 0.017(11) ha 2 

+ 0.017(13) ha 2 \n(ha 2 ) - 0.0035(75) (ha 2 ) 2 \n(ha 2 ) + 0((ha 2 ) 2 ) (24) 
A = 5.42(6), 

for k, = 100. Here the errors are larger because of the greater irregularity among the 
40 random lattices. The goodness of the fit is manifested by the fact that the next- 
to-leading coefficients, (c{ and c\ for i > 1) are clearly decreasing for both k, = 1.3 
and k = 100. The 16 values of Ao(ha 2 ) together with the results of the fit, Eqs. (23) 
and (24), are shown in figure 3 for both values of k. We point out again the manifest 
difference between the two cases. 

However, notice that the coefficients proportional to the divergence in Eqs. (23) 
and (24) are compatible within errors. This is welcome because this coefficient is pro- 
portional to the first terms of the beta and gamma functions, (3$ and 70. For a regular 
square lattice this coefficient is [16,17,23] — l/(47r) = —0.07957, in good agreement with 
the values showed in Eqs. (23) and (24). This is the numerical result which supports 
the fact that the two-dimensional O(N) nonlinear <r-model regularized on a random 
lattice shares the same universality class than regularized on a regular square lattice. 

Now, using Eqs. (20), (23) and (24) into Eq. (21), one obtains the bare Green 
function. We can renormalize it in the MS-scheme 



T^ MS (p 2 ,g r ,h r ^) = Z n TlH e (C,g,h,a)\ , (25) 

h = h r Z g Z~ 1 l 2 
i = a p 

where ri 2 ^ MS is 

r ^^ + (p2 + ^ M _L ln( £, <*» 
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Writing the renormalization constants as Z g = 1 + z g g and Z w = 1 + z w g, we obtain 
z g = (N- 2)(^-ln(£) - c - c' \n(h r a 2 )) - I, 



1 ^ > (27) 

z w = (iV - 1) ( — ln( tj-) - c - c' \n(h r a 2 ) 



Now, using the definition (3 L = —a dj da g = —flog 2 — flig 3 — ■ ■ ■ for the beta function 
and j L = a d/ da In = jog + jig 2 + . . . for the gamma function, we get 

A, = -2c' {N - 2), 70 = -2c (iV - 1). (28) 

Finally, using well known methods [24,25], we obtain the ratio between renormaliza- 
tion group invariant mass parameters, Arandom/ ^-regular- To perform this computation, 
we used the exact values of Co = 51n2/(47r), c' = — l/(47r) and A = 4 for a regular 
square lattice and we assumed the same value of c' for random lattices as well in order 
to obtain finite results for the ratio. This ratio is 



K = 1.3 ^random/ ^regular = exp (-0.042(18) - ~^-^) 

K = 100. ^random/ ^regular = expf 0.19(7) + 0-41 ( 1 ) 



(29) 



N -2 

For the 0(3) non-linear <r-model, these expressions become 0.92(2) and 1.8(2) respec- 
tively. These numbers differ each other more than four standard deviations. Had we 
used the definition Ai, then similar values would have been obtained: 0.92(2) and 2.1(2) 
respectively. The ratios we have got can be compared with those of ref. [26] obtained 
with simulations which give ^-random/ -^-regular ~ 1 .29 5 (in this reference the correlation 
length is measured on rather small lattice and no finite size analysis is done). 

The fact that by varying the k, parameter, one obtains different regularizations of 
the same theory opens some prospects to improve the efficiency of simulations. They 
will be examined in the next section. 



6. Discussion and conclusions 

A family of random lattices has been proposed as a regularization scheme for a 
field theory. This family is characterized by a parameter k, that gives a measure of 
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the "randomness". We have checked that the correct continuum limit is reached for 
a two-dimensional O(N) non-linear <r-model for two quite different values of k. The 
correctness of this statement is measured by the coefficients of the logarithm in Eqs. 
(23) and (24) which should be equal to —0.07957. Indeed, the first universal terms of 
the (3 and 7 functions are proportional to this coefficient. 

We have also evaluated the ratio between lambda parameters for a random lattice 
and a regular square lattice. We have discovered that this ratio depends on k. This 
result means that lattices with different k, are different regularizations. Hence, if an 
average among several random lattices is to be done, one must choose a fixed value 
of k, for all of them. On the other hand, this dependence on k, can also have useful 
consequences. Let us begin by discussing them when Arandom > ^-regular- In this case, 
the scaling window shifts towards larger values of the bare coupling g. Recalling that 
the roughening transition is absent on random lattices (for theories where it has no 
physical meaning) [6,9], one can imagine a scenario where strong coupling expansions 
become useful in this scaling window. If single observables are analysed, all conclusions 
are subject to the fulfilment of the asymptotic scaling condition. 

Let us analyse the case A regu i ar > Arandom- If this happens, then the scaling 
window is placed in a region of lower values of g. This location could favour the onset of 
asymptotic scaling, although a rigorous numerical study is needed in each case to control 
the non-universal terms in the (3 function. Working with lower coupling constant may 
also simplify the subtraction of perturbative expansions from the Monte Carlo signal. 

Actually, the ratio A. r a,ndom/ ^-regular for the O(N) non-linear <r-model in two di- 
mensions has proven to be quite close to 1. Thus, the previous considerations are of 
rather marginal interest for this model. However, in QCD they could be much more 
relevant. Indeed, previous simulations [4,5] have shown that the ratio of lambda para- 
meters is very small: C^IO" 1 ) and O(10~ 2 ) for SU(2) and 577(3) QCD respectively. 
Hence, the K-dependence of this ratio can be more dramatic. In particular, for these 
values of the ratio, the scaling window is strongly shifted towards the region of small g, 
opening the second scenario outlined above. 

The fact that the ratio of lambda parameters is much less than 1 in QCD can be 
easily understood by geometrical arguments [27]. The size of a lattice is usually given 
in units of sites. However, this has no much sense in QCD because for a gauge theory, 
the fields live on the links. So, one should give the size of a lattice as the number of 
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links. On a random lattice in four dimensions there are four times more links than on 
the corresponding regular square lattice with the same number of sites. This explains 
why the lattice spacing on a random lattice is longer than on a regular square lattice. 
For a scalar theory, the degrees of freedom live on the sites and the previous argument 
does not work. 

A numerical study of the consequences of the K-dependence in QCD is in progress. 
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Figure Captions 

Figure 1. Two random lattices with 10 2 sites. The lattice a) was obtained with 
k, = 1.3 and the b) with k, = 100. For both of them we used the random number 
generator RANF and the triangularization method of ref . [3,9]. 

Figure 2. Feynman diagrams contributing to the 1 loop propagator. The cross 
stands for the measure vertex. 

Figure 3. Ao versus ha 2 for k = 1.3 and k = 100. The crosses represent the 
numerical results of table 1. The solid and dashed lines are the curves derived from the 
fits, Eqs. (23) and (24) respectively. 



Table Caption 

Table 1. Values of Ao(ha 2 ) for k = 1.3 and 100. 
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Table 1 



ha 2 


A (ha 2 ,K = 1.3) 


A (ha 2 ,K = 100) 


0.10 


0.46006(4) 


0.4301(1) 


0.15 


0.42573(4) 


0.3950(1) 


0.20 


0.40103(4) 


0.3700(1) 


0.25 


0.38165(4) 


0.3504(1) 


0.30 


0.36566(4) 


0.3344(1) 


0.35 


0.35203(4) 


0.3208(1) 


0.40 


0.34015(4) 


0.3091(1) 


0.45 


0.32960(4) 


0.2987(1) 


0.50 


0.32012(3) 


0.2895(1) 


0.55 


0.31151(3) 


0.2811(1) 


0.60 


0.30362(3) 


0.2735(1) 


0.65 


0.29634(3) 


0.2665(1) 


0.70 


0.28958(3) 


0.2600(1) 


0.75 


0.28328(3) 


0.2540(1) 


0.80 


0.27738(3) 


0.2484(1) 


0.85 


0.27183(3) 


0.2432(1) 
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Figure 1 




Figure 2 




a) 
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